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Quasiparticle dynamics on the topological surface state of Bi 2 Se3, Bi 2 Te3, and superconducting 
CurBi2Se3 are studied by 7 eV laser-based angle resolved photoemission spectroscopy. We find 
strong mode-couplings in the Dirac-cone surface states at energies of ^ 3 and ^15-20 meV, which 
leads to an exceptionally large coupling constant A of ~ 3, which is one of the strongest ever reported 
for any material. This result is compatible with the recent observation of a strong Kohn anomaly 
in the surface phonon dispersion of Bi2Se3, but it appears that the theoretically proposed "spin- 
plasmon" excitations realized in helical metals are also playing an important role. Intriguingly, the 
~ 3 meV mode coupling is found to be enhanced in the superconducting state of CuccBi2Se3. 
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Topological insulators (TIs) are a new class of mate- 
rials with Dirac fermions appearing on the surface pQ. 
The nature of Dirac fermions has already been actively 
studied in the graphitic materials [2 , and it has been 
elucidated that the Dirac band dispersion is anomalously 
renormalized by such effects as electron-phonon interac- 
tion, electron-hole pair generation, and electron-plasmon 
coupling, leading to various intriguing properties [3HS]. 
While topological insulators are essentially understood 
within the noninteracting topological theory [6j [7], the 
Dirac fermions realized in real materials would be af- 
fected by nontrivial many body interactions, and hence 
the investigation of the quasiparticle dynamics is impor- 
tant for extending our understanding beyond the nonin- 
teracting regime. Since the Dirac fermions in TIs are dis- 
tinct form of those in graphitic materials in terms of their 
helical spin texture as well as possible interactions with a 
separate bulk electronic state, the low-energy excitations 
in the topological surface state are of particular inter- 
est. Indeed, such excitations are important not only for 
understanding many body interactions and couplings to 
other degrees of freedom in the topological surface sate, 
but also for assessing the stability of putative Majorana 
fermions that are expected to emerge on the surface of 
superconducting TIs [8] [9] . 

In this context, there are already indications of the 
significance of many body interactions in the topolog- 
ical surface state. For example, a pronounced Kohn 
anomaly to indicate a strong electron-phonon coupling 
was recently observed in the surface phonon branch of 
Bi2Se3 [10]; scanning tunneling spectroscopy (STS) un- 
covered an intriguing feature with finely-resolved sharp 
peaks at low energies (< 20 meV) in the Landau-level 
spectra [TT], pointing to an anomalous increase in the 
quasiparticle lifetime near the Fermi energy (Ep). The- 
oretically, it has been proposed that a novel low-energy 



collective mode called "spin-plasmon" would emerge as a 
consequence of the spin-momentum locking in the topo- 
logical surface state [12]. Therefore, it is important to 
elucidate how the Dirac dispersion is renormalized close 
to Ej?. However, so far the angle-resolved photoemission 
spectroscopy (ARPES) has not been able to detect any 
significant renormalization in the Dirac dispersions in TIs 
[T3] , possibly because of the lack of sufficient energy res- 
olutions. 

In this Letter, we demonstrate that the Dirac disper- 
sion of the topological surface state is indeed anomalously 
renormalized, by using state-of-the-art ARPES with a 7- 
eV laser photon source. Availability of the ultra- high 
energy resolution (~ 1 meV) and the extremely low tem- 
perature (~ 1 K) [14] enabled us to detect low-energy 
kinks in the dispersion at ~3 and ~ 15-20 meV, giving 
evidence for hitherto-undetected mode couplings. The 
analysis of the kinks leads to the estimate of the coupling 
constant A of as large as ~ 3, which is one of the largest 
reported for any material [15j [16] . Despite the existence 
of this boson-mode coupling, we observed no overall band 
reconstruction down to the lowest temperature, indicat- 
ing that the topological surface state is protected from 
density- wave formations, which is usually expected to oc- 
cur with such an extremely strong coupling with bosons 

hzh2T]. 

Single crystals of Bi2Se3 and Bi2Te3 were grown by 
melting stoichiometric amounts of elemental shots in 
sealed evacuated quartz glass tubes. Superconducting 
samples of Cuo.24Bi2Se3 with T c of 3.5 K and a shielding 
fraction of ~ 30% [see Fig. 3(a)] were prepared by elec- 
trochemically intercalating Cu into the pristine Bi2Se3 
[22H25] . ARPES measurements were performed using a 
Scienta R4000 hemispherical analyzer with an ultravio- 
let laser (hv = 6.994 eV) at the Institute for Solid State 
Physics (ISSP), University of Tokyo [26]|27]. 
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FIG. 1: (Color online) Data for (al)-(a5) Bi 2 Se 3 and (bl)-(b5) 
Bi 2 Te 3 . (al), (bl) MDC at E F . (a2), (b2) Band dispersion 
map along T-K [dashed lines in (a3) and (b3)]. (a3), (b3) 
Fermi surface map. (a4), (b4) MDC-derived band dispersion. 
The same dispersion over a wide energy range is shown in 
the lower- left inset. The upper-right inset plots the energy 
difference from the linear dispersion. (a5), (b5) MDC peak 
width Ak (full width at half maximum, FWHM), and the 
Ak(E) multiplied by the group velocity vg(E). The inset 
shows the Ak(E) close to Ep. 

Figures l(al)-(a5) and l(bl)-(b5) show the ARPES 
data of Bi2Se3 and Bi2Te3, respectively. The Dirac cones 
are clearly seen in the dispersion maps [Figs. I(a2) and 
l(b2)], and the shapes of the Fermi surface (FS) are very 
different between the two compounds [Figs. I(a3) and 
l(b3)]. In the present experiment, we did not observe any 
quantum-well states which emerge when adsorption of 
residual gases on the sample surface causes charge doping 
[28H34] . Also, our data are free from spectral intensity 
from the bulk conduction band, as can be clearly seen in 
Figs. l(al) and l(bl) where the momentum distribution 
curves (MDCs) at £V show only two sharp peaks from 
the Dirac dispersion. 

The novel feature in our data is that the MDC-derived 
band dispersions [Figs. I(a4) and l(b4)] obviously devi- 
ate from straight lines, pointing to a large mass enhance- 
ment; the renormalized slope of the dispersion close to 
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FIG. 2: (Color online) Data for Bi2Se3 within a narrow range 
near Ef marked in Fig. I(a2). (a) ARPES image; parabolic 
bands (dashed lines) with a mass of 0.83ra e and 0.14ra e [35] 
are superimposed, (b) Typical EDC with three features, (c) 
EDCs near Jgf on the occupied-state side, and (d) those di- 
vided by Fermi function at T — 7 K. Energies of shoulder-like 
structures [circles in (b) and bars in (d)] are plotted on (a). 
The energy eigenvalue for each k, e(k), which is determined 
by the energy positions at which the spectral intensity be- 
comes maximum, are indicated in (d) with bold bars, and 
plotted in (a) with thick filled hexagons, (e) EDCs beyond 
&f (unoccupied-state side). The dashed line indicates the en- 
ergy of the spectral dip. 



Ep is shown by dashed lines. As shown in the upper 
insets of Figs. I(a4) and l(b4), we calculate the energy 
difference between the putative linear dispersion (which 
is expected when the mode couplings are absent) and the 
measured one to estimate the strength of the coupling as 
a function of energy. We found anomalies at two energies, 
~ —15 and ~ — 3 meV, indicative of electron couplings 
with two different kinds of collective modes. 

The effects of the couplings should also be observed 
in the energy dependence of the MDC peak width (Afc) 
because of the Kramers-Kronig relation between ReE and 
ImE (E is the self-energy). In Figs. I(a5) and l(b5), we 
plot the obtained spectrum of Ak for Bi 2 Se3 and Bi 2 Te3, 
respectively; as expected, Ak(E) presents kinks at the 
two energy scales, ~ —(15-20) and ~ — 3 meV, which are 
better seen in the insets. The ~ — 3 meV kink marks 
the onset of an anomalous increase in the magnitude of 
Ak(E) toward E F . 

To understand the complex behavior of the band dis- 
persions revealed at low energy, we examine the shapes of 
the energy distribution curves (EDCs) around /cp shown 
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FIG. 3: (Color online) Data for Cuo.24Bi2Se3 superconduc- 
tor with T c = 3.5 K and the pristine Bi2Se3. (a) Field-cooled 
(FC) and zero-field-cooled (ZFC) data of the superconducting 
shielding fraction of the sample used for ARPES experiments. 
Band dispersion map along T-K for (bl) Cuo.24Bi2Se3, (cl) 
aged and (dl) fresh surfaces of Bi2Se3- EDCs of Cuo.24Bi2Se3 
close to /cf measured (b2) above T c and (b3) below T c . The 
dashed line in (b3) indicates the spectral dip. EDCs of the 
pristine Bi2Se3 close to k F for (c2) aged and (d2) fresh sur- 
faces. 



in Fig. 2(c). In Fig. 2(d), those original EDCs are 
divided by the Fermi function at the measured temper- 
ature of 7 K convoluted with the experimental energy 
resolution, to remove the effect of Fermi cut-off. In the 
resulting curves, one can identify up to three shoulder- 
like features [an example for k x — k^P c = 0.018A is 
shown in Fig. 2(b)], and the energy positions of those 
features are plotted on the ARPES image shown in Fig. 
2(a). One can see in Fig. 2(d) that one of the shoulder- 
like features on the curves actually corresponds to the 
maximum; the dispersion of this maximum close to hp is 
also plotted in Fig. 2(a) with thick filled symbols, and 
this dispersion is fitted with a parabolic dashed line in 
Fig. 2(a), giving a significantly enhanced effective mass 
of 0.83m e (ra e is the free-electron mass). For compari- 
son, we also plot a putative band dispersion with a mass 
of 0.14m e , which was estimated for the bulk band from 
quantum oscillations [35 . The mass enhancement real- 
ized in the topological state is remarkable. Even more 
surprisingly, the estimated value of A = vo/v^ — 1 (vo 
and are the bare-electron velocity [36] and the renor- 
malized Fermi velocity, respectively) is as large as ~ 3 as 
demonstrated in Fig. 2 (a), which is one of the strongest 
coupling ever reported in any material [TSJ [16] . 
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FIG. 4: (Color online) Results for Bi2Te3. (a) Scattering rate 
(1/t) estimated from EDC peak width As (full width at half 
maximum, FWHM) as well as from the MDC peak width 
multiplied by the group velocity, Ak • vg- Insets show the 
ARPES map at -150 meV (left) and E F (right), (b) Estima- 
tion of As by fitting double Lorentzians to the symmetrized 
EDCs. (c) EDCs far from &f; the energy eigenvalue for each 
k, s(k), corresponds to the peak position in each EDC. 



Perhaps the most direct way to demonstrate the strong 
mode-coupling is to present the peak-dip-hump structure 
in the EDCs. In Fig. 2(e), we show EDCs at fcp an d 
beyond, where the peak-dip-hump shape is usually em- 
phasized, and indeed, a clear dip can be seen at ~ —16 
meV (dashed line). However, a similar structure associ- 
ated with the ~ 3 meV mode is not clearly visible in the 
data of Bi2Se3. 

Intriguingly, we found that a peak-dip-hump structure 
coming from the ~ 3 meV mode coupling becomes visi- 
ble in the superconducting samples of Cu-doped Bi2Se3. 
Figures 3(b2) and 3(b3) shows the EDCs of Cuo.24Bi 2 Se 3 
with T c = 3.5 K measured above and below T c , respec- 
tively. The peak-dip-hump structure is seen at ~ 3 meV 
below T c but it is gone above T c . We also measured the 
pristine Bi2Se3 at the same condition (T = 1.5 K), but 
did not observe the peak-dip-hump [Fig. 3(d2)]. Further- 
more, we doped the surface of the pristine sample up to a 
doping level similar to that of Cuo.24Bi2Se3 by exposing 
it to residual gases [compare Figs. 3(cl) and (bl)], but 
again, the peak-dip-hump structure is not observed [Fig. 
3(c2)]. Obviously, the enhancement of the ~ 3 meV mode 
coupling has something to do with the superconductiv- 
ity, and the origin of this enhancement needs to be scru- 
tinized in future studies. In passing, we note that we did 
not detect any signature of the superconducting energy 
gap in the present experiment, probably because of the 
nonuniform surface condition after the sample cleaving. 
More elaborate studies would be required to nail down 
the topological superconductivity in Cu^Bi2Se3 [24] by 
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ARPES experiments. 

Now we show the relevance of the mode couplings to 
the scattering rate (1/r). Since the MDC width (Afc) is 
not exactly the same as 1/r for a curved band-dispersion, 
we have extracted the scattering rate with two methods: 
one simply uses the peak width of the EDCs (As = 1/r), 
the raw data of which are shown in Figs. 4(b) and 4(c); 
the other calculates the product of the MDC width and 
the experimentally obtained group velocity (Afc • vq « 
1/r). The behavior of Afc • vq in Bi2Se3 and Bi2Te3 
are plotted in Figs. I(a5) and l(b5), respectively, which 
clearly reflect the mode couplings at ~ 20 and ~ 3 meV. 
Intriguingly, as shown in Fig. 4(a) for Bi2Te3, both A^ 
and Afc • vq present a sharp reduction toward E^. It is 
useful to note that this behavior is consistent with the 
STS result showing a sharpening of Landau-level peaks 
at low energies (< 20 meV) [TT] . 

One may also notice in Fig. 4(a) that 1/r gradually de- 
creases toward the Dirac point. This is unusual, because 
a monotonic increase in the electron-electron interaction 
with increasing binding energy is usually expected in con- 
ventional metals [37] . We speculate that this unusual be- 
havior is a consequence of the fact that the penetration 
depth of the surface state increases as the momentum 
moves away from the Dirac point [38 , which makes the 
surface state to gradually gain some bulk character. In 
fact, a feature to suggest such a variation is seen in the 
energy-contour maps plotted in the inset of Fig. 4(a): 
the contour near the Dirac point (left image) is almost 
circular, but close to (right image), it exhibits a C3 
modulation which reflects the symmetry of the bulk. 

Now we discuss the most crucial question, namely, the 
origin of the bosons causing the anomalies at ~ — (15- 
20) and ~ — 3 meV in the ARPES spectra. A plausible 
candidate for the higher binding-energy one is the out-of- 
plane optical phonon mode Af g with uo = 21 and 16 meV 
for Bi2Se3 and Bi2Te3, respectively [39j 00]. It seems 
that Afc begins to decrease toward Ep [Figs. I(a5) and 
l(b5)] at almost the same energy as that of the A\ g mode. 
Also, the relevance of the phonon coupling is supported 
by the fact that the mode energy observed in Bi2Se3 (~ 
20 meV) is higher than that in Bi2Te3 (~ 15 meV), which 
is consistent with the mass difference between Se and Te 
(VrriTe/mse = 1-27). 

As for the ~ 3 meV mode, there are two possible ori- 
gins. One is the optical mode of surface phonons. Re- 
cently, a strong Kohn anomaly was detected by a helium 
atom surface scattering (HASS) experiment in a phonon 
branch of Bi2Se3 [10 at approximately 2fcp irac (fcP irac is 
the Fermi momentum on the Dirac cone) with the char- 
acteristic energy of ~3 meV, and this Kohn anomaly was 
attributed to the surface optical phonon mode [10] . How- 
ever, Kohn anomaly is typically accompanied by an ex- 
tensive FS nesting [20j [4TJ [42] , which is absent in Bi2Se3 
[see Fig.l(a3)]. Also, the anomaly observed in the HASS 
experiment was isotropic [10 , while the Kohn anomaly 



is usually very anisotropic [20l STJ [42] . Hence, the Kohn 
anomaly in Bi 2 Se3 is quite unusual. In our ARPES data, 
the ~ 3 meV mode coupling in Bi2Te3 is obviously weaker 
than that in Bi2Se3, even though in Bi2Te3 the FS nest- 
ing on the Dirac cone is more pronounced [21] and hence 
the Kohn anomaly is expected to be stronger. In addi- 
tion, the A value of 0.43 has been obtained in the HASS 
experiment on Bi2Se3 for the relevant phonons [43 , but 
this value is obviously too small to account for the very 
strong coupling observed here for the ~ 3 meV mode. 
Therefore, the surface optical phonons alone are obvi- 
ously not sufficient for understanding the lower energy 
mode, and we need to seek for additional ingredients. 

In this respect, another, more promising, origin of 
the ~ 3 meV mode is the theoretically proposed "spin- 
plasmon" , which is suggested to have a maximum energy 
of ~ 2.2 meV [12]. This mode consists of coupled plas- 
mons and spin waves, and unlike the Kohn anomaly, it 
is expected for the round Fermi surface as in Bi2Se3 [T2] . 
The observed upturn in the MDC width toward Ep [see 
insets of Figs. I(a5) and l(b5)], which is not expected in 
a typical electron-phonon coupling, can be interpreted to 
signify an increasingly stronger interactions of the Dirac 
quasiparticles with spin-plasmons near Ep. Note that 
such strong interactions between the two are expected 
only when the plasmon spectrum does not overlap with 
the continuum of electron- hole excitations [3] [44], and 
hence the plasmon coupling should dominate the scat- 
terings with q ~ 0. Therefore, it is natural that the 
quasiparticle scattering is enhanced toward Ep in this 
spin-plasmon scenario. All told, it is most likely that the 
large angle scattering (q ~ 2fcp) by the surface optical 
phonons and the small angle scattering (q <C 2fc F ) by the 
spin-plasmons both are playing roles in the enormously 
strong mass enhancement observed near Ep on the topo- 
logical surface state. 

In conclusion, we have investigated the quasiparti- 
cle dynamics in the topological surface state of Bi2Se3, 
Bi2Te3, and Co£Bi2Se3. We found strong mode- 
couplings at the binding energy of ~ 15-20 and ~ 3 meV. 
The coupling to the Af g phonons is proposed as the can- 
didate for the former mode. As for the ~ 3 meV mode, 
there are two possible origins. One is the optical mode of 
surface phonons. The other is the spin-plasmons, which 
are theoretically proposed as low-energy excitations of 
the helically spin-polarized Dirac fermions. Intriguingly, 
despite the extremely large mass enhancement factor A 
of ~ 3, the topological surface state remains free from 
any band reconstruction down to the lowest temperature, 
indicating that the helical Dirac cone is protected from 
density- wave formations which are naturally expected for 
a system with extremely strong couplings to bosons. 
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